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ABSTRACT 

I derive a relation, both 'forward' and 'inverse', between the density and the diver- 
gence of the peculiar velocity which results from the Zel'dovich approximation. My 
calculations assume Gaussian initial conditions. The forward relation expresses the 
density (strictly speaking, the expectation value of the continuity density given the 
velocity divergence) in terms of the velocity divergence, while the inverse relation ex- 
presses the velocity divergence in terms of the density. The predicted scatter in the 
relations is small, hence the inverse relation is close to, though not identical with, a 
mathematical inversion of the forward one. The forward relation is equivalent to the 
well-known 'standard' density-velocity relation in the Zel'dovich approximation. The 
inverse relation, however, is successfully derived for the first time and constitutes a 
potentially interesting alternative to an inverse relation derived by Chodorowski et 
al., based on third-order perturbation theory. Specifically, it may better recover the 
peculiar velocity from the associated density field, when smoothed over scales as small 
as a few megaparsecs. 

Key words: galaxies: clusters: general - galaxies: formation - cosmology: theory - 
large-scale structure of Universe. 



1 INTRODUCTION 

In the gravitational instability scenario for the formation of 
structure in the Universe, the peculiar motions of galaxies 
are tightly related to the large-scale mass distribution. The 
comparison between the density and the velocity fields can 
serve as a test of the gravitational instability hypothesis and 
as a method for estimating the cosmological parameter f2 
(Dekel et al. 1993). In linear regime, the relation between 
the density and the velocity fields is 



CL rigorously computed the mean 5(x) given "&{x), i.e., 
(5)\d, up to third order in (Eulerian) perturbation theory 
(hereafter PT), assuming Gaussian initial conditions. The 
resulting formula is 



2 \ 



(3) 



5(x) = -/(f2) _1 V -v(x), 



(1) 



where f(Q) ~ f2 6 and I express distances in units of 
kms -1 . This equation is applicable only when the density 
fluctuations are small compared to unity. However, sampling 
of galaxies in current redshift surveys and random errors in 
peculiar velocity catalogs enable reliable dynamical analysis 
with smoothing scale of several hT 1 Mpc, where fluctuations 
slightly exceed the regime of applicability of linear theory. 

Relation (|l|) has been recently extended for the mildly 
non- linear regime by Chodorowski & Lokas (1997; hereafter 
CL). Let us define a variable proportional to the velocity 
divergence, 



-/(nr x v • v(x) . 



(2) 



where e| is the variance of the field The coefficients di en- 
tering the above expansion are given by some combinations 
of the joint moments of S and $ and were explicitly calcu- 
lated by CL. Gaussian initial conditions are also assumed in 
the present paper. 

Mildly non-linear relation between the density and the 
velocity divergence is, in contrast to linear relation (|l]), non- 
local. The local estimator of density ^ has thus a non-zero 
variance. Therefore, to obtain an unbiased inverse estimator, 
i.e., of the velocity divergence from the density, we cannot 
simply invert expression (|^). The inverse estimator was ex- 
plicitly constructed up to third order in PT by Chodorowski 
et al. (1998a; hereafter CLPN), who also computed the ex- 
pected scatter in the relation. 

Having said that, it may seem difficult to understand 
why to investigate the density versus velocity-divergence re- 
lation (hereafter DVDR) in the Zel'dovich approximation 
(Zel'dovich 1970; hereafter ZA). This approximation is first 
order in Lagrangian PT and therefore provides only par- 
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tial answers for higher-order perturbative contributions to 
the density contrast and the velocity divergence. Having 
solved the problem rigorously, why to resort to approximate 
schemes again? 

There are a few reasons for which the ZA is still worth 
studying. Firstly, due to its simplicity, it is very popular 
and in wide use. In particular, the density-velocity relation, 
resulting from an Eulerian version of this approximation 
(Nusser et al. 1991), is used in the potent reconstruction 
of the mass density from peculiar velocity data (Sigad et 
al. 1998). (Strictly speaking, Sigad et al. use a formula based 
on the ZA, with the coefficients slightly adjusted to best fit 
N-body results.) It is therefore interesting to see how the 
ZA-based estimator of density relates to perturbative for- 
mula (|). 

Second, N-body simulations have shown that the ZA is 
apparently quite successful in recovering the density from 
the corresponding velocity field. An estimator of density 
resulting from the ZA (the continuity density) happens 
to be merely slightly biased, even for smoothing scales as 
small as a few Mpc (Mancinelli et al. 1994, Ganon 
et al. 1998), where PT is expected to break down. For a 
Gaussian smoothing length of 5 /i -1 Mpc, the rms fluctu- 
ation of a density field is already close to unity and in 
the perturbative expansion for S and terms of all orders 
become comparable. Indeed, for smoothing scales smaller 
than 5 /i _1 Mpc, the value of the coefficient a 2 estimated 
from N-body starts to deviate significantly from the pre- 
dicted value (Chodorowski & Stompor 1998). This is not a 
problem for formula (^), which is applicable to the mass 
density reconstruction from peculiar velocities, a part of 
so-called density-density comparisons. These comparisons 
(e.g., IRAS-potent) currently employ Gaussian smoothing 
length of 12 h -1 Mpc. For such a smoothing length, one 
may hope formula Q to be even better estimator of density 
than the formula based on the ZA. Velocity-velocity com- 
parisons, however, employ smoothing lengths as small as 5, 
or even 3, h' 1 Mpc (e.g., Willick & Strauss 1998). For such 
small scales, an inverted version (i.e., an estimator of veloc- 
ity from density) of the ZA may prove to do better than 
a formula based on rigorous third-order PT. It is therefore 
important to invert the ZA and to test its performance, rel- 
ative to a third-order perturbative formula, against N-body 
simulations. 

Finally, any galaxy density field is derived from a red- 
shift survey, i.e. given originally in the redshift space. To 
compare the galaxy density field with the real-space mass 
density field inferred from peculiar velocity data, the galaxy 
field must be first reconstructed in the real space. The red- 
shifts of galaxies differ from the true distances by the pe- 
culiar velocities, induced themselves by the fluctuations in 
the density field. Hence, the real-space galaxy density recon- 
struction requires a self-consistent solution for the real space 
density and velocity fields. The velocity field remains irrota- 
tional when smoothed over large enough scales, so given the 
field 1}, defined in equation (0), and appropriate boundary 
conditions, the velocity is 

v(x) = ^P- [ d 3 x'0(x') . X ' - * . (4) 
An J \x' - x\ A 

To proceed further, we need a local estimator of the veloc- 
ity divergence from density. Thus, even in density-density 



comparisons, an inverse estimator is indispensable. In linear 
regime # = <5, hence 

V ( x ) = ^ / d 3 x'5(x') , X ' ~ . (5) 

K ' A-K J V '\ X ' -X\* K ' 

Yahil et al. (1991) and Strauss et al. (1992) describe an iter- 
ative technique of simultaneously solving for the real space 
density and velocity fields, in which they use equation (^J). 
However, in the present version of the IRAS-potent com- 
parison, Sigad et al. (1998) include nonlinear corrections 
to this equation. The nonlinear formula for the velocity di- 
vergence in terms of the density they use is a purely phe- 
nomenological fit to CDM N-body simulations. If the ZA is 
applied to predict density from velocity (the potent recon- 
struction), why not to apply it to the inverse case as well, 
i.e. to predict velocity from density (the IRAS reconstruc- 
tion)? The reason why Sigad et al. do not do this is simply 
that thus far, nobody has succeeded in inverting the ZA. For 
example, Nusser et al. (1991) tried to invert it, but failed. 

In the present paper I express density in the ZA as a 
function of the velocity scalars: the expansion (divergence) 
and the shear. This enables me to derive easily the 'forward' 
DVDR in the ZA, i.e., an analog of formula (^). This also 
helps me to invert the ZA, i.e., to compute the mean velocity 
divergence given the density contrast in the ZA ('inverse' 
DVDR). Such an estimator of the velocity divergence from 
density has a scatter, but I show the scatter to be inevitable 
if the estimator is to be local. Moreover, I explicitly compute 
the scatter and find it to be small. In a follow-up paper, we 
test both perturbative and derived here, ZA-based DVDRs 
against N-body simulations (Chodorowski et al. 1998b). 

The paper is organized as follows: in Section ^ I ex- 
press the density in the ZA as a local function of the veloc- 
ity scalars. In Section ^| I average this expression to obtain 
the mean density given the velocity divergence. In Section || 
I compute the expected scatter in this forward DVDR. In 
Section [B| I derive an inverse DVDR, i.e., the mean velocity 
divergence directly in terms of the density. Summary and 
conclusions are given in Section []. 



2 DENSITY IN TERMS OF THE VELOCITY 
SCALARS 

In a Lagrangian approach to PT (Moutarde et al. 1991, 
Bouchet et al. 1992, Bouchet et al. 1995), instead of ex- 
panding the density contrast, one expands the trajectory of 
a particle, 

x = q + DiP (1) (q)+D 2 ^ 2) (q) + . . . . (6) 

Here, q is particle's unperturbed Lagrangian coordinate, x is 
its final (Eulerian) position, D(t) is the linear growth-factor 
of density fluctuations, and ip^(q) and V' (<?) are the cor- 
responding values of the displacement fields ip^ and tp^ 2 K 
The point of the ZA is that only the field V is retained. 
Then, the velocity of a particle at an Eulerian position x is 
simply 

v = f(Cl)H^ 1 \q), (7) 

where H is the Hubble constant, and, expressing distances 
in units of km s _1 , in Eulerian space equation (p) takes the 
form 
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q(x) = x - f 1 v(x) 



(8) 



Hence, the continuity equation reads (Nusser et al. 1991) 

5{x) = \\dq/dx\\ - 1 = ||I- f'^v/dxW -1, (9) 

where the double vertical bars denote the determinant and 
I is the unit matrix. Expanding the determinant in powers 
of products of velocity derivatives we have (following the 
notation of Sigad et al. 1998) 

six) = -r x v • v + r 2 a 2 + r 3 A 3 , 



where 

3 

A 2 (x) = (viji 



A 3 (x) 



and Vi 



i<3 
3 

E 



(10) 



(11) 



(12) 



: dvi /da 



In this paper, by the 'density in the ZA' I always mean 
the continuity density. It is different from the dynamical 
density in the ZA (i.e., resulting from the equation of mo- 
tion), since the ZA conserves mass and momentum simul- 
taneously only to first order. The dynamical density, Sd, is 
exactly as in linear theory, 8a = — / _1 V ■ v. Moscardini et 
al. (1996) used the dynamical density (specifically, the re- 
sulting solution for the velocity in terms of the ZA-predicted 
density) to model the velocity field of clusters of galaxies. A 
large radius, 20 Mpc, of a Gaussian window with which 
they smoothed the velocity field makes this approximation 
indeed applicable. Density-velocity comparisons, however, 
employ considerably smaller smoothing lengths, where the 
linear DVDR is no longer valid. As already stated, one ap- 
proach to find a mildly nonlinear extension of the linear 
DVDR is to rigorously derive higher-order perturbative cor- 
rections. A complementary, less rigorous, but more intuitive 
and also promising approach is offered by the ZA, since the 
continuity density in the ZA successfully recovers the true 
density in N-body simulations by Mancinelli et al. (1994) 
and Ganon et al. (1998). 

Before shell crossing, the velocity field remains irrota- 
tional. This implies that the velocity deformation tensor is 
symmetric and we can decompose it into expansion, 9, and 
shear (the traceless part), aif. 



Vi.j = \98ij + aij , 
where in general 
= § (vij + Vj,i) 

and 



V • v 



Vk,k ■ 



(13) 



(14) 



(15) 



Here, the symbol 5y denotes the Kronecker's delta. Note 
that 9 = —/'!?, "d being defined by equation (0). I will now 
use decomposition ([l3]) in expressions (jl^)-([12]), using the 
methods and the results of Chodorowski (1997; hereafter 
C97). For irrotational fields, the quantity A2 equals to the 
quantity m„ introduced by Gramann (1993). (Note that an 
expression for A2 in Sigad et al. 1998 has the wrong sign.) 



C97 showed that m„ = 6> 2 /3 - 
vorticity term equal to zero), 



- 0-72 (eq.[37] of C97 with the 
where a 2 is the shear scalar, 



a 2 — OijOij (16) 
and I use Einstein's summation convention. Hence, we have 

a 2 = M* 2 -§- 2 )- ( 17 ) 

Substitution of decomposition ([l3|) in equation ([li]) yields 



,1/2.2,2 \ 
+ 3 (°"12 + °"13 + °"23 — 0-11(722 — 0-HCT33 — (T220-33 J 

1 a 3 

where I have used the property an — 0. By definition, 
a 2 = 2 (a 2 2 + <7 2 3 + 0-33) + On + a 2 2 + 0-33 . 
Using the above equation and the identity 



(18) 



(19) 



(an + (T22 + 0-33) =0, (20) 

we can cast the second term in equation ( |li| ) to the form 
a 2 9/6. We thus obtain 

5 (x)=-r^ + |/- 2 (9 2 -ia 2 ) 

-3 I II 



+ r 



(21) 



or, using the variables ■& and 
Sij = — / aij , 
<5(z)=tf + i(tf 2 -§£ 2 ) + ||2 
where 



(22) 

|£ 2 #+itf 3 , (23) 



(24) 



Thus, the density contrast is a function of three scalars, con- 
structed from the velocity deformation tensor: the expansion 
scalar (the velocity divergence), the shear scalar, and the 
determinant of the shear matrix. The above equation is our 
starting point to derive the DVDR within the ZA. 



3 DENSITY IN TERMS OF THE VELOCITY 
DIVERGENCE 

The ZA yields expression ( |lo| ) for the density in terms of the 
velocity derivatives. The resulting expression for the mean 
density in terms of the velocity divergence obtains by aver- 
aging both sides of equation ([Hi]) given the velocity diver- 
gence. Having transformed this equation to the form (p3|), 
the conditional averaging is straightforward. We have 

(S)U =#+l(# 2 - |(S 2 >|,) + (llE^H)U - !<£ 2 >M 



~ 27 17 



(25) 

The Fourier transform of a shear component is (£ij) fe = 
{kikj — \o~ij) 9k, where ki = ki/k and 9^ is the Fourier 
transform of the velocity divergence field. Hence, 



(2tt) 3 



kikj — -Sij ) P#(fc) 







(26) 



(P^(fc) is the power spectrum of the velocity divergence 
field). This means that the shear components are uncorre- 
cted with the velocity divergence. In the case of Gaussian 
random variables, and only in this case, it is a sufficient 
condition to be statistically independent. Since the initial 
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conditions are assumed to be Gaussian, in linear regime 
$ and Eij are independent. When the fields become non- 
linear, they become non-Gaussian as well (e.g., Bernardeau 
et al. 1995, Lokas et al. 1995). In the ZA, however, since 
the velocity field is proportional to the initial displacement 
field, equation (Q), it remains linear even when the density 
field becomes non-linear. Derivatives of a Gaussian field, be- 
ing its linear combinations, are also Gaussian, so •& and Ey 
remain Gaussian, thus independent. In effect, we can simply 
replace the conditional averages in equation j25| ) by the or- 
dinary averages. We have <E 2 ) = (2/3)4- The mean value 
of the determinant of the shear matrix is zero. Hence, 



(*>l* = 


= a[ ZA) d + a[ ZA) (tf 2 -4)+af A) tf 3 , 


(27) 


where 








a[ ZA) 


= 1 - 




(28) 


a[ ZA) 


i 

— 3 ' 




(29) 


and 








a (ZA> 

a 3 


1 

— 27 " 




(30) 



The mean density given the velocity divergence in the 
ZA is thus a third order polynomial in the velocity diver- 
gence, similarly to the third-order PT result (Q). Also the 
coefficients of the polynomial are in many aspects similar to 
the corresponding coefficients resulting from perturbative 
calculations: they form a hierarchy <i ZA) <C <C a Z , 

they are independent of f2, and a~ ZA] has a corrective term,Q 
which scales linearly with the variance of the velocity diver- 
gence field. (The corrective term is due to the term — |S 2 i9 
in equation E3J, a third-order mixed term in the shear and 
the velocity divergence.) Quantitatively, however, the coef- 
ficients are different. As stated earlier, in a separate paper 
we use N-body simulations to test relative accuracy of both 
approximations (Chodorowski et al. 1998b). 



4 SCATTER IN THE RELATION 



Expression ( |27[ ) for density in terms of the velocity di- 
vergence, since obtained by conditional averaging of equa- 
tion (|23j), has clearly a scatter. The rms value of the 
scatter at a given value of the velocity divergence, s|#, 
is given by the square root of the conditional variance, 

\2\|V2 



a|*= ((5-{S)M 2 )\ J . We have 

((s-(s)M 2 )l= + Pdl) 2 ) 



where 



(E 2 >. 



(31) 



(32) 



For 'typical' fluctuations, the first term in parentheses in 
equation ( fjl] ) is of the order of 4 , while the second and the 
third are 0{e%). In large N-body simulations, however, one 
can trace statistical events of the velocity field which are 



* The value of a~ ZA ^ is at the leading order correctly unity, 
because in the limit of small fluctuations the ZA recovers linear 
theory. 



many standard deviations away from the mean. In particu- 
lar, Chodorowski et al. (1998b) reliably estimate the scatter 
as a function of $ even for d well above unity. Therefore, in 
equation (^) I do not assume $ to be small. Since $ and E 
are statistically independent, we obtain 



((S - (S)Uf) \ t = \{y 2 ) (l + \$) 2 + 0(4) . 



(33) 



I recall that given Gaussian initial conditions, the velocity 
field in the ZA remains Gaussian. For such a field, 



(y 2 ) = <(s 2 



(see C97 for details), hence 
2 



2 



2 4 

46 £ <> 



(l+ftf) +0(4)- 



(34) 



(35) 



The probability distribution function for the velocity di- 
vergence has an abrupt cutoff at ■& = —1.5, as PT pre- 
dicts (Bernardeau 1994) and N-body simulations confirm 
(Bernardeau & van de Weygaert 1996, Chodorowski & 
Stompor 1998). Therefore, 1 + ~<d is always positive and 
we obtain finally 



* 



4 , , (ZA) 2 a 



where 



and 



jiBA) _ 1 
°1 _ 9 



0.21 



0.07. 



(36) 



(37) 



(38) 



Carrying calculations up to second order in PT, C97 
derived a formula for a scatter in the DVDR similar to the 
first term in equation (^). Extending the calculations up 
to third order, CLPN derived a formula already containing 
the second term, but were unable to predict the value of the 
coefficient 61. The ZA predicts, in a simple way, not only 
formula ( |3rj ) , but the values of both coefficients bo and 61 as 
well. 

The rms value of the scatter relative to the rms value 
of the divergence, e#, vanishes in the limit — > 0, as ex- 
pected. More importantly, however, this ratio is substan- 
tially smaller than unity even for e$ close to unity. Thus, 
even when almost fully nonlinear, the density and the veloc- 
ity divergence at a given point remain strongly correlated. 
In a follow-up paper we test the prediction of the ZA for a 
scatter in the DVDR against N-body simulations. 



5 VELOCITY FROM DENSITY 

Equation ( jj^ ) can be perturbatively inverted to express the 
velocity divergence as a local function of the density and the 
shear. The resulting expansion for # has an infinite number 
of terms. Up to cubic terms, it is 



i(5 2 -|E 2 )-||E lJ ||-iE 2 5+A 5 3 + 0(e 4 )i (39) 



where e 2 = (<5 2 ). Obviously, the velocity divergence is not a 
function of the density alone. Thus, like the forward relation 
studied in Section (|3j), a local estimator of the velocity di- 
vergence from the density will inevitably have a scatter. To 
obtain an expression for the divergence exclusively in terms 
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of the density, I will average the above equation given the 
density. We have 

= S - ±<5 2 + ±(E 2 )|, - mA\)u - ^ 2 )U5 + ^ 

+0(4) • (40) 

Unlike ■& and E, 8 and E are not independent, since the 
evolved density is a (mildly) non-Gaussian variable. The 
calculation of (S 2 ) |s is a non-trivial problem; I present it 
in Appendix A. The result is 



( e2 )I<5 = §4 



,e 2 s 6 + 0(e 4 s ). 



(41) 



Thus, the mean value of the shear scalar given the density 
weakly depends on the density. The term generating the de- 
pendence is of the order of e|, higher than the constant term. 
This is because at the linear order 8 and E are independent. 
For the same reason, (||£y||)|« = (||Ey||) + 0(4) = 0{ef). 
Usin g th is fact and substituting equation ( |4l| ) in equa- 
tion (E) yields 



where 

JZA) _ 



„(za) 



and 

JZA) 



r[ ZA h + ri ZA) (S 2 - 4) + 4 ZA} 8 3 + 0(4) , (42) 



J-/ 2 

45=-<S i 



(43) 
(44) 

(45) 



And if we instead invert formula (|27|), which expresses 
the mean value of density directly in terms of the velocity 
divergence? Straightforward inversion of ( p7j ) yields expan- 
sion (^), with the coefficients n^ ZA ^ which I will call the 
naive ones, 



JZA) 



AZA) 



n 2 2 

2a 2 £s 



12 
9 t<s • 



(46) 
(47) 



and 



(ZA) i n 2 5 

n\ = -a :i + 2a 2 = . 



The coefficients n. 



(ZA) 



and r^ ZA ^ are equal to r. 



r\ ZA \ respectively, but n\ 



(ZA) 



(48) 
( 2 ZA ^ and 

(ZA) 



different from r[ . This 
is a consequence of a scatter in the DVDR. If a relation 
between two random variables has a scatter, in general the 
inverse relation is not given by a mathematical inversion of 
the forward relation (e.g., forward and inverse Tully- Fisher 
relations). CLPN showed the true and the naive coefficients 
to be related in the following way: 



n = m + (26§ - b 2 )el , 



V2 = n 2 



and 



r-j, = n 3 



(49) 
(50) 



(51) 



Here, bo and 62 are the coefficients entering the leading- order 
perturbative formula for the rms value of the scatter in the 
DVDR, 



s|« = &o4 [1 + 6^7(^4)] 1/2 + 0(4) 



(CLPN; note a slightly different notation used here). This 
expression was derived under an assumption that $ is of the 
order of e#, so the second term under square root is of the 
order of unity. The formula does not account for the second 
term in expression ( prj[ ) because it is already of third order 
in Comparing expressions (^2|) and ( pS6| ) we find that in 
the ZA the coefficient 

b (ZA) = (53) 

and b\ ZA ^ is given by equation (^). From equations 

and 



([so|), ( |5l| ) and (^3]) we have that indeed r. 



(ZA) 



(ZA) 



(ZA) 



Using equations ([£§), ©, @ and (||) 
1 — j^e 2 , in agreement with equation (^ 



obtain 

Thus, I rederived the values of the coefficients r^ A ^ in a 
different way. 

Both 'forward' and 'inverse' relations, expressions ( |27| ) 
and describe mean statistical properties of the matter 
field. They were derived by constrained averaging over all 
possible realizations of the density and velocity fields. It is 
instructive to compare them to the results obtained assum- 
ing spherical symmetry of perturbations. In this case, all the 
shear terms in equation (^3|) vanish, and it simplifies to the 
form 



8 = (1 + i9/3f -1- 
This equation is easily invertible, 
= 3[(l + <5) 1/3 -l] . 



(54) 



(55) 



(52) 



The above expression is in agreement with a result of 
Bouchet et al. (1995) for a spherical top-hat (eq. [A31] of 
Bouchct et al. 1995). When expanded, it yields the values 
of the coefficients r 2 and rz equal to these given by equa- 
tions ( [44] ) and (ffoj). It does not, however, predict the cor- 
rection to the leading-order value of the linear coefficient r\. 
More importantly, it does not involve a constant term. Such 
a term, — r^^e 2 , is present in equation ( |42[ ) and naturally 
assures the (ordinary) mean of the velocity divergence to 
vanish, the property which expression (^i|) lacks. The shear 
terms are thus generally important. 



6 SUMMARY 

I have derived the mildly nonlinear DVDR as predicted by 
the ZA, as well as a scatter in it. The 'forward' relation 
states that the mean density contrast, given the velocity 
divergence, is a third-order polynomial in the velocity di- 
vergence. This is 'a law of Nature' in the ZA, or, rather, 'a 
law of the ZA'. In contrast, the 'inverse' relation, express- 
ing the mean velocity divergence in terms of the density 
contrast, has an infinite number of terms; I have explicitly 
computed the coefficients of the first three. A relation be- 
tween two mildly nonlinear variables should be described 
by a polynomial of third degree quite well. In any case, I 
will not attempt to calculate higher-order coefficients before 
testing the already computed ones against N-body simula- 
tions: the ZA is only an approximation, and modelling its 
prediction even more accurately is no guarantee of bringing 
us any closer to the truth. 

The ^-dependence of the DVDR in the ZA enters only 
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via a factor /(fl), used in the definition (Q) of the scaled ve- 
locity divergence. Similarly, in PT, the relation between the 
density and the scaled velocity divergence is practically f2- 
independent (Bernardeau 1992, Gramann 1993, CL, CLPN; 
cf. also Nusser & Colberg 1998). 

I have explicitly computed a scatter in the forward re- 
lation; a scatter in the inverse relation can be computed 
analogously. I have computed only the 'forward' scatter be- 
cause here we are mostly interested in the mean relations 
and the scatter is only an auxiliary quantity informing us 
about the limitations of our local estimators. The predicted 
scatter is relatively small, even for the fields which are al- 
most fully nonlinear. Therefore, the inverse relation, when 
obtained by a straightforward inversion of the forward, will 
be only slightly biased. Indeed, a proper calculation of the 
inverse relation yields a minor correction to the value of the 
linear coefficient and no corrections to the quadratic and 
cubic coefficients. This offers an efficient way of deriving 
approximate values of the coefficients of the higher-order 
terms, if there is in future any need to include them. 

I have not included the effects of smoothing the evolved 
density and velocity fields, while the fields inferred from ob- 
servations are smoothed. In rigorous PT, smoothing slightly 
changes the values of the 'forward' and 'inverse' coefficients, 
making them weakly dependent on the underlying power 
spectrum of mass fluctuations and the window function used 
(CL, CLPN, Chodorowski et al. 1998b). The inclusion of 
smoothing in the ZA can be done in an analogous way to 
that in PT. Formula (^o|), which is clearly for unsmoothed 
fields, is however found a successful estimator of smoothed 
density from smoothed velocity in N-body simulations by 
Mancinelli et al. (1994) and Canon et al. (1998). Thus, the 
apparent success of the ZA is somewhat accidental: neither 
higher orders in Lagrangian PT nor the smoothing effects 
are included, but it still works. If it really works, an in- 
verse estimator based on the ZA may, on scales smaller than 
about 5 ft _1 Mpc, do better than the corresponding esti- 
mator based on third-order (Eulerian) PT. If so, it would 
be very useful in large-scale velocity-velocity and density- 
density comparisons. In a follow-up paper, we will test the 
performance of both approximations in N-body simulations 
(Chodorowski et al. 1998b). 

The perturbative derivation of the mildly nonlinear 
DVDR by CL and CLPN is very formal. In contrast, a local 
relation between the density and the velocity scalars in the 
ZA (though strictly valid for unsmoothed fields only) en- 
ables one to derive the DVDR in an easy and intuitive way. 
In this picture, the source of the scatter in the DVDR is the 
second, 'hidden', parameter of the velocity field, namely the 
shear. (In reality, another source of the scatter is smooth- 
ing.) A similar conclusion was drawn by C97 on the basis of 
second-order PT. However, in the ZA-formula for the den- 
sity there are terms up to third order in the velocity deriva- 
tives. The third-order terms give rise to the second term of 
formula ( |36| ) for the scatter, unpredictable by second-order 
formalism. Furthermore, a third-order, mixed term in the 
shear and the velocity divergence is a source of the correc- 
tion to the leading-order value, unity, of the linear coeffi- 
cient, ai, in the forward relation. The conditional average of 
this term, given the velocity divergence, yields an additional 
term linear in the velocity divergence, with the coefficient 
proportional to the variance of the velocity divergence field. 



The linear scaling of the correction to ai with the variance 
is indeed formally predicted by PT. Thus, even if the ZA 
fails to model quantitatively the mildly nonlinear DVDR, it 
will remain a useful tool to understand it. 
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APPENDIX A: CONDITIONAL AVERAGE OF 
THE SHEAR SCALAR 

I outline here a derivation of the expectation value of the 
velocity shear scalar given the velocity divergence or the 
density contrast. The only assumption made is that the den- 
sity and the velocity fields remain in the mildly non-linear 
regime, so the calculation can be performed perturbatively. 
Besides that the derived formula is entirely general, i.e. it 
can be applied to any approximation of mildly nonlinear 
dynamics. Here, I apply it to the ZA. 

The derivation of (E 2 )|,j is greatly simplified when we 
introduce an auxiliary variable 



= + E 2 - {E 2 } , 



(Al) 



where <E 2 ) is an ordinary average of the shear scalar. Ex- 
panding shear components in a perturbative series, Ey = 
E« + £< 2) + yields E 2 - <E 2 > = E (1)2 - <E« 2 } + 
2E 1 < 1) E I (2) + 0(4), where E^ 2 = E^E^. Similarly, = 
0(i) + tf( 2 ) + 0(») + 0(e%), or, for short, = 0i + 2 + $3 + 
0(4)- Hence, 



= fix + fa + As + e>(4) ■ 



(A2) 



In the expressions above, e is the linear variance of the 



velocity divergence field, e 



i), and the symbol 



stands for the connec ted (reduced) part of the moments. 
From equations (A6) and (A7) we have 



<E 2 )|„ = (£ 2 } + ( Cl -l)tf + c 2 (0 2 



- 4) + c 3 tf 3 • 



(A14) 



Using expansion (A2) of in the expression for the coeffi- 
cient C2 yields 



2e 4 c 2 



(t»iGfti-0a)) 



0« 2 (e (1) 



0. 



(A15) 



In the last but one step I used the fact that and E are 
independent at linear order (see Section ^) . Casting similarly 
the coefficients ci and C3 we obtain 



<E 2 )|„ = (E 2 ) + si40 + s 3 3 + 0(4) , 
where 



(A16) 



where 

01 = di , 

02 = 2 + E 2 , 

and 

03 = tf 3 + E 3 , 



-(i)" 



2E (1) E (2) 



By definition (Al), 
(E 2 )|^<E 2 ) + </3)U-0, 



(A3) 
(A4) 

(A5) 
(A6) 



so our problem reduces to the calculation of mean given 
0, where both variables, of vanishing mean, are mildly non- 
linear and equal to each other at linear order. This problem 
was solved by CL. According to CL, 



{0)U = dtf + C2 (tf 2 - 4) + c 3 tf 3 + 0(4) , 

where 

(S3/3 — SW) Sstf Z4 



Cl 



C'2 



C3 



1 + 



Z 2 + ■ 



2 



S30 — i?3fl 

6 ' 
Z4 — (S33 — S 3 #) S30 
6 



(A7) 

(A8) 

(A9) 
(A10) 



(I use here slightly different notation.) The quantity S3/3 is 
defined by 



e*S 3 p = H0102) ■ 



(All) 



and is defined in an analogous way. The quantities _Z 2 
and Z± are given by 



e"Z 2 = {0 2 $2) - (n) + (0103) 
and 



3> 



e 6 Z 4 = 3(/3 2 /3 2 2 ) c - 3(0 2 0^) c + {0i0 3 )c 



(A12) 
(A13) 



Sl = Z2 — ~Z± ! 



S3 



6 



with 



£ Z2 



and 



(l? 2 E 2 ) + (l?iE 3 



e 6 ^ 4 = 3(0 2 2 E 2 ) c + (0?E 3 ) c 



(A17) 
(A18) 

(A19) 
(A20) 



Thus, the average value of the shear scalar given the velocity 
divergence is equal to its unconstrained average plus the 
corrective terms, dependent on the divergence. These terms 
are of the order of 4> higher than the first term, which is of 
the order of e|. This is again due to the fact that at linear 
order, E and are statistically independent. 

Though the variable denotes the velocity divergence, 
the only specific property of I have used thus far was its in- 
dependency of E at linear order. This property is also shared 
by the variable 5, since Si = $1. Therefore, an expression for 
the mean value of the shear scalar given the density contrast 
can immediatel y be w ritten by replacing the symbol with 
5 in expression dAlq). Specifically, 



(E 2 ) \ s = (E 2 ) + s' lS 2 S S + s'zS 3 + 0(4) , 
where 

s i ~ Z 2 — -Z 4 , 

' - 1 cr' 

s 3 — a ^4 ' 

6 

with 

4,Z 2 = (&Ea) + <5iE 3 ) 
and 

e 6 ^ = 3(5 1 2 52E 2 ) c + ((5 1 E 3 > c . 



(A21) 

(A22) 
(A23) 

(A24) 
(A25) 
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The formulas (A16) and (A21) are general in a sense that 
they are applicable to any approximation of mildly nonlinear 
dynamics (including rigorous PT). Here, I will apply them 
to the ZA. 

In the ZA, the velocity field remains linear all the time, 
so i?2 ZA ' = Sjj 2 ^ = 0. (The quantity Ej 2 ^ is non-zero 
because it is constructed from first-order quantities) . Hence, 

Z (ZA) = Z (ZA) = Q 

(E 2 )|f A) = (E 2 ). (A26) 

This result is otherwise obvious, since in the ZA, E and i9 
remain independent in the nonlinear regime (see Section^). 

The case of (E 2 ) |^ ZA ' , however, is not so trivial, because 
unlike the velocity field, the density field in the ZA is non- 
linear. From equation (|23|), 

4 ZA) = | (V 1)2 - |S (1)2 ) ■ (A27) 
This yields 



.£_l var ( E (D 2 ) , (A28) 

where var(E« 2 ) is the variance of the linear shear scalar. 
This variance is given by equation ([m]), hence 

Z^ ZA) = -±. (A29) 
45 

Analogous calculation shows that 

Z' A (ZA) = (A30) 

(I recall that Z± is constructed from the connected part of 
moments). This yields s'^ 2 ^ = -4/45 and s' 3 (ZA) = 0, 
hence 

(E 2 >|f A) = (E 2 )-| e 2 5 + 0( e l). (A31) 

The ordinary average of the shear scalar is equal to (2/3)e| 
and the variance of the velocity divergence field is not equal 
to the variance of the nonlinear density field, Eg. The differ- 
ence, however, is C(e|). Therefore finally 

(E a )|f A) = |d-^ + 0(4). (A32) 
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